"Splitting" of Delocalized States in a Double— Layer System in a Strong Magnetic Field 
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A double-layer system in a strong perpendicular magnetic field is considered. We assume a 
random potential in each layer to be smooth. We also assume that there is no correlation between 
random potentials in different layers. Under these conditions the equipotential lines from different 
layers, corresponding to the same energy, may cross each other. We show that, if the tunnel coupling 
between the layers exceeds some characteristic value (which is much smaller than the width of the 
Landau level), then the probability for an electron to switch equipotential (and, thus, the layer) at 
the intersection is close to one. As a result, the structure of each delocalized state in a double- 
layer system becomes completely different from that for an isolated layer. The state is composed 
of alternating pieces of equipotentials from different planes. These combined equipotentials form a 
percolation network. We demonstrate that the regions, where equipotentials from different planes 
touch each other, play the role of saddle points for such a network. The energy separation between 
two delocalized states is of the order of the width of the Landau level, and the critical exponent of 
the localization length is 7/3-the same as for an isolated layer. 



PACS Numbers: 72. 20. My, 73.20.Jc, 73.40. Hm 



I. INTRODUCTION 

In a recent paper |l| , S0rensen and MacDonald studied 
the integer quantum Hall effect in a double-layer system 
with uncorrelated disorder. They investigated numeri- 
cally the localization properties of single-electron states 
for two limiting cases: of a short-range and of a smooth 
disorder. For a short-range disorder, the relevant param- 
eter is t/Y, which is the ratio of the tunnel integral and 
the disorder-induced width of the Landau level. If this 
ratio is large, the situation is transparent since localiza- 
tion of states, belonging to symmetric and antisymmetric 
combinations of the size-quantization wave functions, oc- 
curs independently. Thus, one has two delocalized levels, 
located at the two well-separated maxima of the density 
of states, and spaced in energy by 2t. 

When t/Y <C 1, the density of states has only one 
maximum; the usage of the basis of symmetric and anti- 
symmetric combinations is inadequate. However, it was 
demonstrated in the simulation of Ref. Q (see also Ref. 
) that the separation of delocalized states is close to 2t 
even for t/Y as small as 0.1. In other words, the positions 
of delocalized states with uncorrelated disorder appear to 
be the same as in the case of correlated disorder. 

To study the case of a smooth disorder, the authors 
[jlj employed the network model proposed by Chalker 
and Coddington ||, which was generalized to describe 
a double-layer system in a way similar to that in Refs. 
||-|| . In the Chalker-Coddington model, the derealiza- 
tion results from the tunneling of an electron, moving 
along equipotential lines, through the saddle points of 
the random potential. The randomness of the potential 
is included by assuming the phase acquired by an elec- 
tron traversing a link (equipotential line) to be random. 
For a double-layer system each link of a network carries 
two channels. The authors tried several variants of incor- 
porating the inter-layer tunneling through the coupling 



between the channels. The results obtained were essen- 
tially the same. In contrast to the case of two spin-split 
levels with random mixing the authors observed 

the doubling of the correlation length exponent. One of 
the main conclusions made in Ref. [Q on the basis of the 
simulation performed, is that there is no observable split- 
ting in the positions of the delocalized states within the 
range of t/Y studied, which makes doubtful the applica- 
bility of the network model to the double-layer systems. 
Motivated by this observation, in the present paper we 
take a microscopic approach to the problem. 

In a strong magnetic field the structure of electronic 
states in the presence of a smooth potential is determined 
by its topography. The crucial difference between the 
cases of correlated and uncorrelated disorder is that in 
the latter case the equipotentials from the different lay- 
ers, may cross each other after projection on the same 
plane. This is illustrated in Fig. 1. As a result of such 
a crossing, there is a finite probability for an electron to 
change equipotential (and, correspondingly, the plane). 
Obviously, this probability increases with increasing t. 
We will show that the probability V, for an electron to 
stay within the same equipotential after crossing, is given 
by 



V = exp 



2nt 2 l 2 



|vi x V2|7i 



(1) 



where / is the magnetic length and vi, V2 are the drift 
velocities of an electron along the first and the second 
equipotentials respectively. These velocities are related 
to the values £\, £2 of the electric field in each of the 
planes at the point of crossing as: V1.2 = c(£\ t 2 xB)/B 2 , 
where B is the magnetic field perpendicular to the planes. 
Then we can rewrite (M as V = 

cxp(— 27rf 2 /e 2 |£i x £2!' )■ This sets a relevant scale for t. 
Indeed, the fields can be estimated as Y/eR c , where 
R c is the correlation radius of a smooth potential. 
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FIG. 1. Crossing of equipotentials from the first 
(long-dashed lines) and the second (dashed lines) layer after 
projection on the same plane. 

Thus we see, that the probability to retain equipo- 
tential is close to 1 if t <C Tl/R c . In the opposite case, 
t ^> ri/R c , practically each crossing results in the change 
of equipotential. It is very important that the crossover 
value t ~ Tl/ R c is much smaller than T if the disorder is 
smooth. This implies that, as t gradually increases, the 
regime when the equipotential is switched with probabil- 
ity close to 1 establishes already at rather small values of 
t, where the density of states still represents a single peak. 
The applicability of the formula (|l|) to the double-layer 
system with uncorrelated random potential is limited by 
the condition t<T. This is because in the derivation of 
(|l|) it was assumed that electric fields £\, £2 are constant 
within the entire region where the "interaction" between 
equipotentials occurs. One can see from the derivation 
presented below, that this region is of the order of the 
magnetic length when t <C eSx^l, i.e. when V is close 
to 1. However, in the opposite limit, t 3> e£i^l, which 
corresponds to the strong coupling, the extension of the 
interaction region is of the order of t/e£i t 2 and is much 
larger than On the other hand, for Eq.(Q) to be rel- 
evant, this extension should be much smaller than R c , 
which leads us to the condition t <C T. Note that for 
t r, the language of symmetric and antisymmetric 
states becomes adequate. Then, similarly to the case of 
the short-range disorder, the density of states represents 
two peaks, each having a delocalized state in the center. 

Thus, we have established that in a wide region 
Tl/R c <C t <C r, each crossing of equipotential lines from 
different planes leads to the change of the plane, in which 
the electron moves. As a result, the saddle points || of 
the random potentials V\{p) and V^(p) in the planes, 



which played a crucial role for derealization at t = 0, 
become irrelevant in this region (a typical saddle point 
would be bypassed due to switching of equipotentials, see 
Fig. 2). 
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FIG. 2. Bypassing of a saddle point due to switching of 
equipotentials at the intersections. 

It appears that within the region Tl/R c < t < V, 
the delocalized states have completely different structure 
than at t = 0. Their energy positions, measured from the 
center of Landau level, are E = ±E C , where E c is of the 
order of the width of the Landau level T (see Fig. 3), and 
they are composed of alternating pieces of compact loops, 
which are equipotentials V\ (p) = E c and V2 (p) = E c for 
one state and Vi(p) = —E c and V2(p) = —E c for the 
other state. To justify this scenario, assume first that 
energy E is deep in the tail of the Landau level. Then 
the equipotentials V\{p) = E and V^(p) = E, being pro- 
jected on the same plane, represent a set of isolated cir- 
cles. As E moves up, the equipotentials, corresponding to 
different planes, start to overlap and form clusters, as it 
is shown in Fig. 4 . It is important that, due to switching 
of equipotentials at intersections, the motion of an elec- 
tron within a cluster occurs either inside the cluster or 
along its boundary (see Fig. 4). As E further increases, 
the average size (correlation radius) of the clusters grows 
and, at some E — E c ~ T, the classical percolation oc- 
curs. As the energy sweeps through E c , critical clusters 
first come close, touch each other, and, finally, merge. 

To describe the corresponding change in the electronic 
states, quantum mechanics becomes important. We will 
show that there is a complete correspondence between a 
region of touching of equipotentials and a conventional 
saddle point. Namely, the transmission coefficient, de- 
fined as a probability for an electron to change the clus- 
ter while passing through the region of touching, has the 
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same form as the transmission coefficient of a saddle point 
I 



where E c is the "height" of the effective saddle point 
and it differs from the value of potential at the point 
of touching by an amount of the order of t. For the 
characteristic energy A we derive below an estimate: 
A ~ T 3 /' 2 l 2 /t 1 / 2 R 2 . It is easy to see that in the region 
Tl/R c <t<T one has A <C I\ 

Thus we see, that the description of derealization at 
E ±E C reduces to a single-channel network model with 
critical clusters, playing the role of links, and the regions 
of touching acting as nodes. 



plane for the electron motion. We can present only a 
plausible argument about the evolution of the positions 
of delocalized states in this limit. Note that, although 
a typical crossing is not efficient for t < Tl/R c , the re- 
gions of touching of equipotentials from different planes 
still act as saddle points. It will be shown below that the 
corresponding condition on t is t > T(l/R c ) 4 / 3 . Thus, 
in the domain Tl/R c > t > r(l/R c ) 4 / 3 , these effective 
saddle points would couple electronic states belonging to 
closed equipotentials from different planes. As a result, 
the delocalized states would occur at energies ±E^\ at 
which the size of a closed equipotential, R C {T / 'E^) 4 / 3 , is 
big enough to have ~ 1 effective saddle point somewhere 
on its perimeter. It is obvious that E^ <C E c . Then the 
problem again reduces to a single-channel network model 
with the cells of the network being closed equipotentials, 
alternatingly, from the first and from the second layer, 
and the nodes being the effective saddle points. 
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FIG. 3. Schematic plot of the energy positions of delocal- 
ized states in a double-layer system at different values of in- 
terlayer coupling. 

This implies that the delocalized states have the same 
critical exponent of the localization length p,p|-^2|, k = 
7/3, as in an isolated layer. The magnitude of the local- 
ization length can be estimated as follows. Quantum me- 
chanical description of the transmission through the ef- 
fective saddle point becomes important when \E — E c \ < 
A. The size of the "unit cell" of the network for E — 
E c ± A is ~ i? c (r/A) 4 / 3 , where 4/3 is the critical expo- 
nent for the classical percolation. Then we have, 



E — E r 



I 2 



"n — 
I e 



R c tf{E-E c )i 



(3) 



As t becomes smaller than Tl/R c , the above picture of 
electronic states is not valid anymore. A typical cross- 
ing of equipotentials would not cause the change of the 




FIG. 4. A cluster formed by equipotentials from different 
layers after projection on the same plane. Full lines show the 
resulting trajectories of the electron motion. 

To estimate the magnitude of E^ note that the 
perimeter of critical equipotential scales with energy as 
{□j: C(E) ~ i? c (r/£:) 7 / 3 . The probability that two 
equipotentials from different layers would come close and 
form an effective saddle point with height within the 
interval (E — A, E + A) can be estimated as A/T ~ 
(r/t) 1 / 2 (Z/i? c ) 2 . Remember that A is the energy scale of 
the transmission coefficient in Eq. (Q). Then the condi- 
tion to find EP can be written as: (£(^ 1) )/E c )(A/r) ~ 



I, and it yields E { c 1} - T{Y /tf/ 1A {l/ R c ) 

In the next two sections we derive the basic equations, 
(|l|) and (||). Section IV concludes the paper. 
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II. TRANSITIONS BETWEEN INTERSECTING 
EQUIPOTENTIALS 

Suppose that two equipotentials from different planes, 
corresponding to the energy E, measured from the center 
of the Landau level, intersect at p = 0. The potentials 
V\{p), V 2 (p) behave near p = as 



V 1 (p)=E + e£ 1 - p; V 2 (p) = E + e£ 2 ■ p 



(4) 



Then the amplitudes \&i and ^> 2 , to find an electron re- 
spectively in the first and in the second plane, satisfy the 
following system of equations 



m \ 



eS\ ■ p — hui n 



+ e£ 2 ■ p — Tito n 
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2/ 


















I 




2y 







'2 



(5) 
(0) 



We chose a symmetric gauge A = B(—y, x, 0)/2 and per- 
form a transformation to the new variables, defined as 



l(X - a), 



... d d 

v = - tl{ dx + dS 



d_ 

dx 



1 

21 



d 

dX 



d_ 
ds~ 



d_ 

dy 



21 



(X + s). 



(7) 
(8) 



This transformations was previously used by Fertig and 
Halperin || to describe the transmission through a sad- 
dle point. It allows to separate effectively the cyclotron 
motion (s-coordinate), and the motion of the guiding 
center (X-coordinate). With new variables, the system 
(H),(0) reads 



ds 
Pi x X 



92 -R 8 



ds 
d_ 

dX 



2n 



PlxS *i 



2t 

1)*! = ^ 



(9) 



d 2 d \ 

-Q^-lp2yQ- s + S 2 -f32 xS )y2 + 



dX 



(10) 



where we have introduced the notations: 0i,2x = 
2e£\^ x l/hjjj, Pi,2y = 2e£i t 2yl/?w and u> — eB/mc, is 
the cyclotron frequency. The term proportional to d/ds 
can be eliminated by the following substitution 



*i, 



exp 



-1 — — s 



01 



$1,: 



(11) 



As a result the equations 
form 



lOh acquire a separable 



<T>, [ PixX - ip\y£x ) l|> '- 



ds 2 ' V '2 

(2n + 1 + f3j) $1 = 

|L exp {i ^-^ a - fahv-^fov ) $ 2 , (12) 



d& 



2 



$2 + IhxX - i& 



■2y 



dX 



(2n + 1 + Pi) $2 = 

2i J .Ply — 02y .Plxflly — 02x1- 

exp < -1 - s + 1 - 



*i, (13) 



where 01 = P 2 x + (3 2 y and /3| = 02x+02y Since X enters 
linearly in both equations, we will search for solutions in 
the following form: 



<!>, = F„[.Vio„ [s-^f 



G n (X)(f) n ( S 



fox 



(14) 
(15) 



where <p n are the oscillator eigenfunctions. Then the sys- 
tem of equations for the coefficients F n , G n , describing 
the motion of the guiding center, takes the form 



-Wiy-Hx + PixX^j F n - f3 2 F n = ^-G n , (16) 
-ifou-Hx + fo**) G « - Pi°n = ^Fn, (17) 



We have set unity the overlap integral of the func- 
tions 4> n centered at 0i x /2 and /?2x/2. Indeed, one has 
Pix - 02x ~ e£ lt2 l/hu> < 1. In fact, the system @,(|TJ) 
is equivalent to the system describing non-adiabatic tran- 
sitions between the crossing energy levels in molecules. 
This problem was first considered more than 60 years 
ago [ [l4j. Up to a phase factor, the solution of the system 
( |l6| ) , (|17|) can be expressed in terms of parabolic cylinder 
functions |TJ, D l/ (±xV2e m / 4 ), with v and X given by 
the following formulas: 



2 72 



t 2 l 



— _ / 2£i v £ 2 y 



e 2 l 2 \£i x £ 2 \ |vi x v 2 |?i : 

hui e n (£i y — £ 2 y) 



X 



2 el\£ix£ 2 \ 



t \£ix£ 2 \ 

Using the asymptotics of the D-functions, 
D V {X) ~X»e-i x \ 



(18) 
(19) 



\X\ 



CO, 



37r „ Sir 

< arqX < — 

4 4 



(20) 



4 



D V (X) ~ X v e~* 
\X\- 



1/2 



(gTT) 

r(-iz) 



00, 



one can see that the right behavior at X — > +00 for the 
function F„(X) (no reflected wave), is insured by the fol- 
lowing choice: 



F n {X) oc D v (±xV2e ln/i 



(22) 



The asymptotics for F at X — > 00 and X — > —00, dif- 
fer by a factor exp(z7ri/). With 1/ given by fllSj), the 
probability to retain the equipotcntial after crossing, 
7> = |F(-oo)| 2 /|F(oo)| 2 , takes the form Eq. (g). 

As it was mentioned in the Introduction, the character- 
istic spatial scale of the interaction region depends of the 
value of t. If t <C eEi^l, we have |z/| <C 1 and the char- 
acteristic X in ( p2| ) is of the order of 1. This means that 
the spatial scale for the interaction is I. In the opposite 
case, t ^> e£\fil, the parameter v is large, \v\ ^> 1, and 
the product X v exp(— X 2 /2) in ( p2| ) sets a scale for X: 
X ~ l^l 1 ^ 2 ~ t/e£i 2 l, which leads to the characteristic 
scale ~ t/e£i } 2 for the interaction region. 



III. THE EFFECTIVE SADDLE POINT 



Assume for concreteness that two equipotentials that 
nearly touch each other, have their origin in two displaced 
potential minima in each layer (Fig. 5) 



TT?0 2 

V 1 (x,y) = ^[(y-y 1 ) 2 +x*]+V 10 



mf2f 



[x 2 - 2y iy + yj] +V W , 



V 2 (x,y) = ^[(y + y 2 ) 2 +x 2 ]+V 2a 



' 2 [x 2 + 2y 2 y + y 2 ] +V20, 



(23) 



(24) 



where (yi,0) and (—y 2 ,0) are the positions of the min- 
ima, Vio,V 2 q are the heights and fii,f2i are the curva- 
tures. We neglect the terms ~ y 2 in (p^), p^ ) since the 
relevant y appears to be small. 



2 

s s 



/ 



\ 



FIG. 5. Effective saddle point formed by equipotentials 
from different layers after projection on the same plane. 

The condition that two equipotentials come close to 
each other at x = 0, y = 0, can be expressed as 

— ^ + V w « — ^ + V 20 ~E (25) 

Then the system of equations for the amplitudes ^i,^> 2 , 
similar to (||), (||), can be written as 



Px 

2m 



m(uj 2 + ttl) 



Py 



uj 2 + Qf m 



( 0? 



V 2mu>'- 



:p y - mCl^iy J *i 



E — Tilo I n + — 



Vu 



*i = W 2 



(26) 



pI_ 

2m 



m(u 2 + 



Py 



(J3 



2 -2 



E-huj 



1 



p 2 + mn 2 y 2 y ) * 2 



-V, 



20 



* 2 =<*i. (27) 



2) 2 

We will search for a solution in the following form: 



= J dkA{k)e lky (j> n (x 
dkB{k)e lky (/) n (x 



Lu 2 + n 2 



uj 2 + q 2 



i 2 k 



l 2 k 



(28) 
(29) 



Upon substituting (g9J) into (|26|), (|27|), we get a 

system of equation for the functions A and _B 
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dA 

- imn\yx— + a(k)A(k) = tc(k)B(k), 
ok 

f) Fi 

+ imflly 2 — + b(k)B(k) = tc(k)A(k), 
ok 

where c(k) is the overlap integral 
c(fe)= 

dxcj). 



(30) 
(31) 



UJ 2 l 2 k 



n 2 



n 2 



and the coefficients a(k) and b{k) are defined as 



a(k) = 


mVLll^k 2 


(- 


mVt\y\ 


2 


2 


b(k) = 


m n 2 2 i A k 2 

2 


(- 


mQ%y% 
2 



- Vi 



U 



20 



(32) 

(33) 
(34) 



ft can be easily seen that c(k) = 1 for f^i = ^2 and 
the correction in the case when the two frequencies are 
different, is proportional to (fif — Vt^) 2 /uj a . Since the 
the random potential is smooth, we can neglect this cor- 
rection and set cik) = 1. The system (|30|), ( |3l| ) can be 
reduced to a single second-order differntial equation, say, 
for A(k) 



d 2 A ( a b 
y. i 

dk 2 \m£l 2 yi mOfjfe 



ab-t 2 



da 



i 2 Q 2 Q 2 ,yiy2 mQ 2 yi dk 



dA 
~dk + 

A = 0. 



(35) 



The term with first derivative can be eliminated by the 
following substitution 



A{k) = exp 



after which Eq. 

d 2 A 
dk 2 4 
t 2 



" dk' ( a ^ 



b(k') 

\m^lfyi mfl 2 y2 



A(k), 
(36) 



']) takes the form 
a b 



mflfyi mfl 2 y2 
il 4 k ( 1 1 

2 \yi V2 



A = 0. 



m 2 Sl\Sl\y\y2 

ft is convenient to introduce the following notations: 
2 _ 1 1 

yo yi z/2 ' 



(37) 



\/yw2 
yo 



(38) 



£0 



2n 2 yi 

^lyo 



mn 2 y 2 ' 

& Kio 



2^2/2 



E 



2 



20 



(39) 



Using the definitions of a and b, the equation (|37j) be- 
comes 



d 2 A 
dk 2 



+ 



l A k 2 



mfigyo 2y 



il 4 k 
Vo . 



.4 = 0. (40) 



Eq. ( |40| ) has the form of the Schrodinger equation with a 
complex "potential energy" . However, not all the terms 
in the "potential energy" are relevant. This becomes ob- 
vious if we introduce the following rescaling of the argu- 
ment 



mQ 2 y 2 



Then the equation (^ 

2 £ . e 



d 2 A 

dz 2 



t 



t 

takes the form 
- 1 2 1 

L , 1 4 2 

-= — V —z a 



IZOt 



(41) 



A = 0, (42) 



where we have introduced the dimensionless parameter 



/yo 



t 



(43) 



This parameter can be estimated as follows. The typi- 
cal value of yo is ~ i? c ; the curvature tto can be found 
from the condition |l6[ ]: mQ^R 2 ~ T. Then we have 
a ~ (//i? c )(r/t) 3 / 4 . But within the domain we are 
interested in, t is larger than Tl/R c . Then we get 
a < (l/Rc) 1 / 4 <C 1. This allows to drop the last two 
terms in the potential energy. It can be also seen that 
the effective saddle point corresponds to |eo + t\ € t. 
Indeed, under this condition Eq. (B2|) takes the form 



d 2 A 

dz 2 



,2 O (eo + *) 
ta 2 



A = Q. 



(44) 



which is the equation, describing the scattering from 
the inverted parabolic potential. The expression for the 
transmition coefficient for this potential is well-known 
E7| 



T(E) 



1 + exp 



27r(e +t) 
to 2 



(45) 



We see that the characteristic energy scale for the change 
of the transmission coefficient is 



rs 



tot r- — 



(46) 



Using the above estimates, the condition a <C 1, which 
guarantees that the region of nearly touching of two 
equipotentials acts as a saddle point, can be rewritten as 
t T(l/ i? c ) 4 / 3 . Note that A decreases with increasing 
t, reflecting the fact that the larger is t, the closer should 
equipotentials approach each other in order to form the 
effective saddle point. 
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IV. CONCLUSION 



The main result of the present paper is that for a 
double-layer system with smooth uncorrelated disorder, 
the position of delocalized states as a function of cou- 
pling between the layers, has a wide plateau (Fig. 3). 
This is in contrast to the case when the disorder in both 
layers is strictly the same. In the latter case, the energy 
separation between the localized states is just 2t. Note 
that the prediction we make can, in principle, be tested 
experimentally on a single sample, since, as it was shown 
in Ref . [Esj , the tunnel integral t can be effectively tuned 
(suppressed) by applying a parallel magnetic field. 

When the authors of Ref. |Q questioned the applicabil- 
ity of the network model to double-layer systems, they 
implied that the network remains the same as the cou- 
pling between layers changes. Our considerations shows 
that a single channel network model applies within a wide 
range of t, but the network itself experiences restructur- 
ing with t. In particular, while at t = the nodes are the 
saddle points of the intralayer potential, at large enough 
coupling the nodes are the points where equipotentials 
from different layers touch each other. 

Experimentally it might be hard to realize both, 
smooth disorder within each layer and the absence of 
correlation between the layers. However, for our argu- 
ments to apply, it is sufficient that the correlation be- 
tween the layers is not absolute. In other words, it should 
be allowed for two equipotentials, corresponding to the 
same energy, to be displaced by more than a magnetic 
length. Then our picture remains valid, but the correla- 
tion would lead to the shift of the percolation threshold 
to an energy much closer to the center of the Landau 
level, than in the case of uncorrelated disorder. This im- 
plies that the "splitting" of the delocalized states in the 
regime when the coupling between equipotentials is al- 
ready strong would still remain much smaller than T, i.e. 
the plateau shown in Fig. 3 would become lower, and 
hence, narrower. 

In the paper we demonstrate that when the tunnel in- 
tegral exceeds the characteristic value Tl/R c , the critical 
exponent, k, of the localization length is the same as 
for t = 0. We also argue that the description based on 
a single-channel network model (and thus k = 7/3) is 
applicable within the interval Tl/R c > t > T(Z/i?c) 4 / 3 . 
Therefore if the doubling of k established in Ref. jj] oc- 
curs, this may happen only for t < T(l/ i? c ) 4 / 3 . However, 
it seems more likely that the description in terms of a 
single-channel network model with neighboring cells be- 
longing to different layers applies even at very small t. 
The reason why we anticipate this is the following. The 
solution of the model problem in Sect. 



Ill shows that 



if two minima of the random potential in different lay- 
ers are located anomalously close to each other (yi , 7/2 in 
(p3|),(|24|)) are anomalously small), then the formation of 
the effective saddle point becomes possible even at very 
small t. Certainly, the smaller is t, more sparse these sad- 



dle points are. At t = the localization length increases 
in each layer as £(E) cx E~ 7 / 3 , and each localized state 
consists of many cells of the intralayer network, separated 
by saddle points. Then, if t is finite and very small, at 

some E = Ec there will be ~ 1 effective saddle point 

per perimeter of a localized state in each layer. At this 

(2) 

energy (as well as at E = —Ec ) the states in two layers 

would form a new network with a much larger unit cell 

(2) 

~ £,{E C '). Then the localization length would behave as 
(E - Ei 2) )- 7 / 3 with a prefactor much bigger than that 
for a single layer. Certainly, this scenario is only hypo- 
thetical. 

Note that there is a significant difference between our 
picture and the one outlined by S0rensen and MacDonald 
0. The line of argument in Ref. [Q is as follows. In the 
absence of a disorder an electron residing initially, say, 
in the first layer, would oscillate between the layers with 
a period r = 2irh/t. When the disorder is present, the 
electron drifts within the layer with velocity v\. It was 
assumed in Ref. [0 that the change of the layers would 
most probably occur after an electron travels the distance 
l<ir = v\T = 2irhvi/t ~ Tl 2 /tR c . In our picture the pos- 
sibility for an electron in the first layer to tunnel depends 
on the actual topography of the random potential in the 
second layer. In contrast to Ref. |Q, in our picture the 
change of the layers occurs locally, at the intersections of 
equipotentials with the same energy. At the same time 
in the regime of strong coupling, t > Tl/R c , that we 
considered, the distance (in the vicinity of an intersec- 
tion), over which the change of the layers takes place, is 
~ tR c /T and it is much larger than ld r - 

Our picture also differs from that by Laikhtman and 
Menashe [l^] . Similarly to Ref. [jj], they assume that 
the process of changing layers, during the drift along 
the equipotential, occurs homogeneously, but they get 
a different estimate for the characteristic length of travel 
within a given layer. Their estimate, R^T /It, is larger 
than tR c /T and, correspondingly, larger than l^ r . 
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